Quantum Hall ferromagnetism in graphene: an SU(4) bosonization approach 
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We study the quantum Hall effect in graphene at filling factors v — and v = ±1, concentrating on 
the quantum Hall ferromagnetic regime, within a non-perturbative bosonization formalism. We start 
by developing a bosonization scheme for electrons with two discrete degrees of freedom (spin-1/2 and 
pseudospin-1/2) restricted to the lowest Landau level. Three distinct phases are considered, namely 
the so-called spin-pseudospin, spin, and pseudospin phases. The first corresponds to a quarter-filled 
[y = — 1) while the others to a half- filled [y = 0) lowest Landau level. In each case, we show that the 
elementary neutral excitations can be treated approximately as a set of n-independent kinds of boson 
excitations. The boson representation of the projected electron density, the spin, pseudospin, and 
mixed spin-pseudospin density operators are derived. We then apply the developed formalism to the 
effective continuous model, which includes SU(4) symmetry breaking terms, recently proposed by 
Alicea and Fisher. For each quantum Hall state, an effective interacting boson model is derived and 
the dispersion relations of the elementary excitations are analytically calculated. We propose that 
the charged excitations (quantum Hall skyrmions) can be described as a coherent state of bosons. 
We calculate the semiclassical limit of the boson model derived from the SU(4) invariant part of the 
original fermionic Hamiltonian and show that it agrees with the results of Arovas and co-workers for 
SU(N) quantum Hall skyrmions. We briefly discuss the influence of the SU(4) symmetry breaking 
terms in the skyrmion energy. 

PACS numbers: 71.10.-w, 81.05.Uw, 73.43.-f, 73.43.Lp 
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I. INTRODUCTION 

Graphene consists of a single atomic layer of carbon ar- 
ranged in a honeycomb lattice. 1,2 When an uniform per- 
pendicular magnetic field is applied, the system displays 
an unconventional integer quantum Hall effect (QHE), 3,4 
where the Hall conductivity a xy — 4(n+ l/2)e 2 //i (n in- 
teger) and the filling factor is defined as v — 4(n + 1/2). 
Such unusual behavior of o~ xy is understood within a 
single-particle model 2 ^ which shows that each Landau 
level in graphene is approximately four-fold degenerate 
(valley, the so-called K and K' points, and electron spin). 

More interesting, experiments performed at higher 
magnetic fields showed new quantum Hall plateaus at 
v = 0, ±1, and ±4, 6 indicating that the degeneracies of 
the n = and n = 1 Landau levels are lifted, fn par- 
ticular, for v — ±4, the behavior of the minimum of the 
longitudinal resistance R xx in terms of the total mag- 
netic field suggests that here the quantum Hall effect is 
due to the lifting of the spin degeneracy of the n = 1 Lan- 
dau levels However, the origin of the plateaus at v = 
and v = ±1 is not completely understood. Different sce- 
narios were proposed. It was suggested that the effect 
is due to Coulomb interaction, which favors a quantum 
Hall ferromagnet ground stated Alicea and Fisher— pro- 
posed that the plateaus might be related to symmetry 
breaking terms, such as Zeeman and underlying lattice 
interactions, which give rise to a paramagnetic phase as 
it occurs at v = ±4. An explanation based on the so- 
called "magnetic catalysis" mechanism was proposed by 
Gusynin et al£ This mechanism predicts that the long- 
range Coulomb interaction generates an excitonic gap, 



which lifts the valley degeneracy only of the lowest Lan- 
dau level. In combination with the Zeeman splitting, the 
observed quantum Hall plateaus at v = and ±1 are un- 
derstood. More recently, Abanin et al. w argued that the 
transport response of the quantum Hall state at v = is 
due to counter-circulating edge states. 

In this paper, we study the quantum Hall ferromag- 
netism in graphene via a non-perturbative bosoniza- 
tion method for the case of electrons with spin-1/2 and 
pseudospin-1/2 restricted to the lowest Landau level. It 
constitutes a generalization of the formalism 11 recently 
proposed by one of us to study the two-dimensional elec- 
tron gas at v = 1 realized in GaAs heterostructuresi 12 i 13 
Within this formalism, the elementary neutral excita- 
tions (magnetic excitons) and the skyrmion-antiskyrmion 
pair excitations of the system are described in the same 
framework, namely an effective interacting boson model. 
Such method is quite general and was used to calculate 
spin excitations of the fractional quantum Hall systems at 
v = 1/3 and v — 1/5^ as well as to study Bose-Einstein 
condensation of magnetic excitons in the bilayer quan- 
tum Hall system at total filling factor vt = 1 (spinless 
case)J£ 

Concerning the latter, the great majority of mod- 
els proposed to study this system assumes fully spin- 
polarized electrons. However, nuclear magnetic reso- 
nance measurement a 16 ' 17 indicate that the electron spin 
degree of freedom might be relevant. Indeed, it was sug- 
gested that the incompressible-compressible phase tran- 
sition observed in this system may involve a modification 
of the spin polarization.™ Therefore, theoretical tools 
which allow us to properly treat the electron-electron in- 
teraction and simultaneously take into account the elec- 
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tron spin and layer (pseudospin) degrees of freedom are 
needed. The formalism developed here might be also use- 
ful to study the bilayer quantum Hall system at vt = 1 
in GaAs heterostructures (spinfull case). 

Our paper is organized as follows. In Sec. II, we define 
the creation and annihilation boson operators and derive 
the boson representation of the (projected) electron den- 
sity, spin, pseudospin, and mixed spin-pseudospin density 
operators. Three distinct cases are considered, the so- 
called spin-pseudospin phase, which occurs when the low- 
est Landau level is quarter-filled, the spin and pseudospin 
phases, which are related to a half-filled lowest Landau 
level. In Sec. Ill, we apply the generalized bosonization 
formalism to study the QHE in graphene at v = and 
v = ±1, focusing on the quantum Hall ferromagnetic 
regime. Our starting point is the effective continuous 
model recently proposed by Alicea and Fisher.— For each 
quantum Hall state, an effective interacting boson model 
is derived and the dispersion relations of the elementary 
neutral excitations arc analytically calculated. We com- 
ment on some possible effects of the boson-boson inter- 
action and show how the quantum Hall skyrmion might 
be described within this scheme. A summary of the main 
results is presented in Sec. IV. 



II. THE BOSONIZATION METHOD 

In order to develop a bosonization scheme for elec- 
trons with two discrete degrees of freedom (spin- 1/2 and 
pseudospin- 1/2) and restricted to the lowest Landau level 
subspace, we follow the lines of Ref. [TT| and start by 
studying the corresponding noninteracting model. 

Let us consider N noninteracting electrons moving 
in the xy plane under a perpendicular magnetic field 
B = Bz. In addition to the electronic spin (a, A =|, J,), 
let us also include a discrete pseudospin index a, (3 = ±. 
Restricting the Hilbert space to the lowest Landau level, 
the kinetic energy is quenched and therefore the Hamil- 
tonian of the system is 

H = H z +Hvz (1) 
= -\H I d " r (^ + oBp)$t ff (r)$ a(7 (r). 

In addition to the Zeeman term 7iz, we also include an 
extra term (TL-pz) which breaks the pseudospin degener- 
acy. As we will see below, a finite Ep helps us to define 
a set of different reference states. Ez = gfJ-BB is the 
Zeeman energy, where g is the effective electron ^-factor 
and hb is the Bohr magneton (see Appendix [XJ. CT (r) 
is a fermion field operator that can be expanded in the 
(Schrodinger) lowest Landau level basis \n = Om) (sym- 
metric gauge as 

*L(r) = ^(n = 0m|r)4 QCT , 



* Q(J (r) = ^2(r\n = Om)c maa . (2) 

tn 

The operator cj noc(7 (c m aa) creates (destroys) an elec- 
tron in the lowest Landau level, with guiding center m, 
pseudospin a, and spin a. Substituting Eq. ([2]) into Eq. 

one sees that the Hamiltonian H is diagonal in the 
lowest Landau level basis, i.e., 

W= -- 53 ^2( aE z + aEp)cl ia(T c maa . (3) 

rn— a,(T 

The above Hamiltonian has four highly degenerate en- 
ergy levels, whose energies are — (Ez + Ep)/2, —(Ez — 
E P )/2, (E z — E P )/2 and (E Z + E P ), and the degeneracy 
of each level is = l/2irl 2 . Here, I = y/hc/eB is the 
magnetic length and we assume that the total area of the 
system is one. In the following, we will concentrate on 
three distinct configurations of the system: total number 
of electrons N = and Ez > Ep, which we call spin- 
pseudospin phase; N = 2N^ and E z > E P (spin phase); 
and N = 2N0 and Ez < Ep (pseudospin phase). 

As discussed in Ref. [ill , the creation and annihilation 
boson operators are defined by considering the neutral 
(particle-hole) excitations above a well-defined reference 
state. As each one of the above phases has a different 
reference state (noninteracting ground state), the three 
cases will be analyzed separately. However, before do- 
ing that, we should firstly discuss the representation and 
the algebra of the electron density, spin, pseudospin and 
mixed spin-pseudospin density operators projected into 
the lowest Landau level. 



A. Density operators and the lowest Landau level 
algebra 

We start by defining the following projected density 
operator 

p a *M*) = *L(r)*/5A(r), (4) 

where the fermion field operators are given by Eq. ^j, 
and whose Fourier transform is 

/WA(q) = J d 2 r e-^ r *L(r)^ A (r) 

= E J d 2 re-^( m \v)(v\m')cl aa c mipx 

m,m' 

= e^ )2/2 GmMW&arCm'PX, (5) 

with q = |q|. The function G mym i(x) is defined as 

n n \ at i \ [™^ f- il (<lx - i( ly) Y" 
Gm,m'(iq) = ( m ~ m ^y^n\ ^ ) 
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xL£'- m ((lq) 2 /2) 



a, i\ m '- ( -H(<h + iQy) 



xi"'((/<z) 2 /2), 



(6) 



where L™r m (x) is the generalized Laguerre 
polynomial. 18 Due to the fact that the operators 
Paa-.j3\(q) are projected into the lowest Landau level, 
their commutation relations are modified, i.e, 

[*WA(q) , PaV', /3 'A<(q')]=e q ' q '' 2/2 

where q A k = l 2 (q x k) • z. 

It is convenient to introduce an isospin index I such 
that 

1 = (a, a) = (+,T), (+,|), (- T), (-,1) - 1,2,3,4. 

In this new representation, the commutator J7J) simply 
reads 

b/./(q),P/j(q')] = e^/^Sjj^'VprM + q!) 

-S^je-^'/'prM + q!)]. (8) 

It is also useful to define a four-component spinor (r) 
as 

*t( r ) = (*t_ T ( r ) *^( r ) *t_ T ( r) $t ^ , (9) 

which, in the isospin language, assumes the form 

$t(r) = (*} (r) *t (r) *t (r) *t (r) ) . (10 ) 



The (projected) electron density operator can now be 
written in terms of the spinor (fT0|) as 



/J (r)=l't(r)#(r)=^*t(r)^(r) (H) 
i=i 

and therefore its Fourier transform may be expressed in 
terms of the density operators pu(q) as 

p(q) = [pn(q) + P22(q) + /?33(q) + /Mq)] ■ (12) 

The same can be done for the spin-cr and pseudospin-a 
electron density operators 

PT(q) = pn(q) +P33(q), Pi(q) = P22(q) + P44(q), 

(13) 

p+(q) = pii(q) + /»22(q), p-(q) = ^(q) + ^(q)- 



The definition (|T0|) implies that the structure of 
the spin-pseudospin space is SU(2)ps <& SU(2)spin 
and therefore, the spin, pseudospin, and mixed spin- 
pseudospin density operators (7i = 1) are respectively 
defined as 

S(r) = ^&(r){l 2x2 ®a)*(r), (14) 



P(r) = -*t( r )( fr(8 i ax2 )$( r ) j (15) 



PS(r) = 4 f (r)(<T®<T)$(r), (16) 

Here, 1-2x2 is the two-dimensional unit matrix and a = 
(a x <j y a z ) is a vector whose components are the Pauli 
matrices. Expanding the Fourier transform of the compo- 
nents of S(r) and P(r) in terms of the density operators 
P/./(q), we have 

^(q) = ^ |>u(q) - P22(q) + ^(q) - P44(q)] , 



S+(q) = [5x(q) + i5 Y (q)] - Pl 2 (q) + Mq), (17) 



5"(q) = [5 x (q) - i5 Y (q)] = P2i(q) + P4 3 (q), 



and 



p z(q) = - [pn(q) + P22(q) - P3s(q) - Pu(q)] , 



^ + (q) =Pi3(q)+p 24 (q), (18) 



p (q) = P3i(q) + P42(q)- 

Similar considerations hold for the mixed operators 
PS(r), in particular, we have 

PzS z {q) = - [pu(q) - p 2 2(q) - A33(q) + /»44(q)] • (19) 

This component of PS(q) will be important in the next 
sections. We should mention that the representation 
([T4")) - (fT())) does not correspond to the standard represen- 
tation of the special unitary group SU(4) [see Ref. [lj^ 
for details], but it follows the ideas presented in the Ap- 
pendix A of Ref. [U 

Finally, with the aid of the commutator 1(5)), a long 
but straightforward calculation shows that the density 
operators (fi"2l) . (fT4|) and (fT5|) obey the lowest Landau 
level algebra (the same results have been derived in a 
more general way2i) 

[p(q),p(k)] = 2zsin(qAk/2) e q - k / 2 p(q + k), 
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- iqAq,/ V/(q+q') 



(22) 



FIG. 1: (color online) Schematic representation of the four 
highly degenerate lowest Landau levels when Ez > Ep. The 
state |SPFM) is obtained by completely filling the energy level 
(f). bf, b2 and b3 are the elementary neutral excitations 
which are related to the density operators P2i(q), P3i(q), and 
P4i(q), respectively. 



[/£(q),p(k)] = 2»sin(qAk/2)e q - k / 2 P;(q + k), 
[7£(q),/£(k)] = (i/2)* , 6 siii(qAk/2)e* k /V(q + k) 
+ ie abc cos (q A k/2) e q ' k/2 ^ (q + k), 

[P a (q),S b (k)} = ^sin(qAk/2)e q • k / 2 P a S b (q + k).(20) 

Here, a,b,c = X, Y, Z and e abc is the Levi-Civita tensor^ 
fi = S,P and therefore /f(q) and 1^ (q) stand respec- 
tively for S a (q) and P a (q). Due to the fact that the 
density operators (fT2")l , (fT4")) , and (fl"5)) are projected into 
the lowest Landau level, the algebra (|20|) is different from 
the usual one of the generators of the SU(4) groupi 19 i 20 

B. Spin-pseudospin polarized state 

Let us now study the noninteracting system described 
by the Hamiltonian ([3]), assuming that the total number 
of electrons N = and E z > E P . The four (highly 
degenerate) energy levels are schematically displayed in 
Fig. [TJ In this case, the noninteracting ground state of 
the system is a spin-polarized pseudospin-polarized state, 



|SPFM)= [] <4 + T |0>, 



(21) 



where |0) is the fermion vacuum. Notice that the neu- 
tral (particle-hole) excitations are created by applying 
the density operators P2i(q), P3i(q), and P4i{q) on the 
state |SPFM). 

From Eq. ©, it follows that the commutator between 
each one of the above density operators and its respective 
Hermitian conjugate is 

bi/(q),pn(q')] = e q - q '' 2/2 [ e 4qAq7 Vi(q + q') 



with / = 2,3,4. By expanding the density operators 
pn(q) around the (reference) state ([2T|) . 



p u (q) = (SPFM|p7j(q)|SPFM) + *pjj(q) 



= N05/ 4 5 q .o + <5p//(q), 



(23) 



and neglecting the fluctuations with respect to the aver- 
age value, the commutator (|2"2"|) assumes the form 



[pi/(q),pn(q')] 



<5 q .- q 'N e^') 2 / 2 . 



(24) 



One can see that, although the relations (|22|) do not 
correspond to the usual canonical commutation relation 
between the annihilation and creation boson operators, 
their expectation values in the ground state |SPFM) do. 
In other words, as long as the number of particle-hole ex- 
citations in the system is small, i.e., (pn(q)) » Spn(q), 
the density operators P2i(q), P3i(q), and P4i(q) may be 
approximately considered as boson operators. Moreover, 
by noticing that 

P32(q) = P42(q) = P43(q) ~ o, 

which is related to the fact that the average values of the 
above density operators with respect to the state defined 
by Eq. pTjl vanish, it turns out that the three kinds of 
boson operators are independent. 

Based on the above analysis, we define the following 
set of creation and annihilation boson operators 



&I(q) 
4(q) 
4(q) 



a?P2i(q), 
a g P3i(q), 
agP4i(q), 



b 2 (q) 
& 3 (q) 



a?Pi2(-q), 

a?pi 3 (-q), 

agPu(-q), 



(25) 



with a q — / 4 / v/N^. From now on, we will assume 
that the above operators obey the usual canonical algebra 

p»t(q) J 6)(k)] = [6 i (q),6 i (k)]=0 J 

[6 i (q),6 J t(k)]=^5 q , k . (26) 

Finally, we should mention that the reference state 
|SPFM) is indeed the boson vacuum as one can easily 
show that 6i(q)|SPFM) = 0. 

Once the boson operators are defined, the boson repre- 
sentation of any operator O is determined by examining 
the commutators [0, b|(k)] (i = 1,2,3) and the action 
of O in the reference state |SPFM). For instance, let us 
consider the density operator pn(q). From Eqs.© and 
(|23|) . we have 

[pH(q)At(k)] = -e-^ 2 / 4 e-^ k / 2 6i(q + k), 
with % = 1,2,3. Moreover, 

Pll (q)|SPFM) = Ny5 q , |SPFM). 
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Using the fact that the three kinds of boson operators 
(|2"5)l are independent, the above relations are satisfied if 
the density operator pn(q) is expanded in terms of the 
bosons bi(q) as 

pu(q) = N^o ~ e'^ 2 / 4 e- 4qAk/2 &l(q + k)ix(k). 

k,i 

(27) 

Similarly, it is possible to show that 



P22(q) 


k 


e iqAk/2 6 t( q . 


-k)6 x (k), 


P33(q) 


= e -(i?) ! /4y 

k 


el qAk/2 6 t (q _ 


-k)6 2 (k), (28) 


P44(q) 


k 


el qAk/2 fo t( q _ 


-k)6 3 (k), 



i.e., the expansions of all density operators /0//(q) in 
terms of bosons are quadratic. 

With the aid of the relations (|27|) and (|28|) . one can 
easily write down the boson representation of the elec- 
tron density [Eq. (fT2|) ]. the z-components of the spin 
[Eq. ([Tf| ] and pseudospin [Eq. IT8|) ] density opera- 
tors, and the mixed spin-pseudospin density operator 
PzS z {(l) [Eq. (0], namely 

p(q) = N <5 q , o + 2ze-^ 2 / 4 

x^sin(qAk/2)foi(q + k)6 4 (k), (29) 

i,k 

J£(q) = 5N^ q ,o + ^/f(q,k)6i(q + k)6 i (k), 



(30) 



with /^(q) = 5z(q),Pz(q), and <Sz-Pz(q), and the form 
factors fi{x) are given by 

/f(q,k) = / 3 s (q,k)=- e -^ 2 / 4 cos(qAk/2), 
/ 2 S (q,k) = Je -^ 2 / 4 sin(qAk/2), 




FIG. 2: (color online) Schematic representation of the four 
highly degenerate lowest Landau levels when Ez > Ep. The 
state |SFM) is obtained by completely filling the energy levels 
(1) and (3). bl, b2, b3, and b4 are the elementary neutral 
excitations which are related to the density operators P2i(q), 
P23(q), P«(q), and P43(q) respectively. 



Appendix [B] . Another important point is that some- 
times the Hermiticity requirement is not full-filled. For 
instance, in the spin-pseudospin phase, the expansion of 
P2i(q) in terms of bosons does not correspond to the one 
of pxi{— q). As it was already discussed in Ref. |Tl|, it 
does not constitute a major problem because the boson 
expressions (J2H]), O, (ED) and |B2]), derived within the 
procedure outlined above satisfy the lowest Landau level 
algebra (j2"0j) . 

The asymmetric boson representation found for some 
operators might be related to the fact that the bosoniza- 
tion method explicitly breaks some symmetries. For in- 
stance, in the spin-pseudospin phase, the spin "direc- 
tions" up and down are no longer equivalent because the 
bosons fo;(q) are defined with respect to the reference 
state |SPFM). As a consequence, the bosonic expressions 
of the spin density operators S~(q) = P2i(q)+P43(q) and 
5 + (q) = pi2(q) + P34(q) [see Eq. (fT7|)] are asymmetric. 
We will see later in Sec. IIIDI that the bosonic representa- 
tion of the spin density operators satisfies the condition 
5 + (q) = [S' _ (— q)]t because for the pseudospin phase 
only the pseudospin symmetry is explicitly broken. 



/f(q,k) = le ^) 2 / 4 sin(qAk/2), (31) 

/ 2 p (q,k) = / 3 P (q,k) = - e -^) 2 / 4 cos(qAk/2), 

and 

/f S (qk) = / 2 PS (q,k)-- e -^) 2 / 4 cos(qAk/2), 

/ 3 PS (q,k) = ie -^ 2 / 4 sin(qAk/2). 

In addition to the set of density operators analyzed 
above, the boson representation of the density operators 
P2i(q) and P34(q) and the respective Hermitian conju- 
gates pi2{— q) and Piz{— q) will be useful in the next 
section, where the bosonization scheme will be applied 
to study the QHE in graphene. Sometimes, the expres- 
sions are not so simple as the one presented above [see 



C. Spin phase 

In this phase, Ez > Ep and the total number of elec- 
trons N = 2N0. The ground state of the noninteract- 
ing Hamiltonian ([3]) is a spin-polarized pseudospin-singlct 
state, 

|SFM)= [] cL- T cL +T |0). (32) 

m— 1 

The particle-hole excitations are now created by the den- 
sity operators pai(q), /023(q), P4i(q), and ,043 (q) as it is 
illustrated in Fig. O 

The commutation relations between p/j(q) {I = 1,3 
and J = 2,4) and their respective Hermitian conjugates 
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PJ/(-q) rea d [see Eq. ©] 

[p/j(q),PJ/(q')l = e q - q '' 2 / 2 



3 4qAq ' /2 P//(q + q') 



-e-^'/ipjjiq + q')]. (33) 

Here, the expansions of pii(q), /»22(q), p33(q), and P44(q) 
around the reference state |SFM) are given by 

PlI (q) = (SFM|p J/ (q)|SFM) + 5p Ji (q) 

= N^,(5/,i+(5j,3)5 q ,o + 5pji(q), (34) 



and therefore the commutation relations (|33p reduce to 
(neglecting the density fluctuations <5p//(q)) 



[p/./(q),p.//(q')] ~ tfq-kN^ 



(35) 



Using the same arguments of the previous section, we as- 
sume that P2i(q), P4i(q), P23(q), and P43(q) are approx- 
imately boson operators. Indeed, they are independent 
operators because 

pu(q) = (SFM|pjj(q)|SFM) + *p/j(q) wO, (36) 

for (I, J) = (1,3) and (4,2). 

To sum up, the spin phase is characterized by a set of 
four independent boson operators defined as 

&{(q) ee a g p 2 i(q), 6i(q) ee a 9 pi 2 (-q), 

4(q) = a g/°23(q), 62 (q) = a q p 32 (-q), 

4(q) = ««P4i(q), 63 (q) = «,Pi4(-q), (37) 

&I(q) = ««/°43(q), 64 (q) = a 9 p 3 4(-q), 

with a g = e^ ? ' / 4 /-i/iN^, and obeying the canonical bo- 
son algebra (|26|) . 

The introduction of new boson operators implies that 
the expansions of the density operators p//(q) in terms 
of the bosons are no longer given by Eqs. (|27| and (|28|) . 
Following the same procedure discussed in the previous 
section, it is possible to show that 

Pll(q) = N^ q , - E e-(^ 2 / 4 -^ k / 2 6j(q + k)6,(k), 




FIG. 3: (color online) Schematic representation of the four 
highly degenerate lowest Landau levels when Ez < Ep. The 
state |PFM) is obtained by completely filling the energy levels 
(1) and (2). bl, b2, b3, and b4 are the elementary neutral 
excitations which are related to the density operators p3i(q), 
P32(q), P4i(q), and P42(q) respectively. 



N0<5 q ,o — * 2N^i5 qi o. However, the boson representation 
of the z-componcnts of the spin and pseudospin density 
operators and the mixed operator PzSz{q) are modified, 
i.e., 

5z(q) - N (5 q ,o - e ~ (/9)2/4 

x ^cos(qAk/2)&|(q + k)6 2 (k), (39) 

i,k 

^(q) = E/H^^blCq + ^ft^k), (40) 
i.k 

where £\ = Y*=v 7 z(q) = p z(q) and S z P z (q), and the 



k i=l,3 



P22(q) 
P33(q) 

P44(q) 



£ e-(^) 2 / 4+ ^ Ak / 2 6l(q + k)6i(k), 

k i=l,2 

N^ q , - E e-^ 2 / 4 -^ Ak / 2 6t(q + k)6 l (k), 

k i=2,4 

^ e -(^) 2 /4+,qAk/2 & t (q + k) ^ (k ^ (38) 



k z=3,4 



form factors are given by 

/f(q,k) = -/f(q,k) = ze-^ 2 / 4 sin(qAk/2), 
/ 2 p (q,k) = -/ 3 P (q,k)= e -^) 2 / 4 cos(qAk/2),(41) 

and 

/f S (q,k) = -/f S (q,k) = - e -^) 2 / 4 cos(qAk/2), 
/ 2 PS (q,k) = -J 3 PS (q,k) = - l e-^ 2 / 4 sin(qAk/2). 

Finally, the new bosonic expressions of Pi2(q), P2i(q), 
P34(q), and P43(q) are shown in the Appendix [B] [sec 
Eqs. jl3| and ijBg])]. 

D. Pseudospin phase 

The situation here is quite similar to the one discussed 
in the previous section, because again N = 2N0 but now 
Ez < Ep. As a consequence, the ground state of the 
nonintcracting model |[3]) is a spin-singlet pseudospin- 
polarized state, 



By adding up the four terms above, one can see that 
the electron density operator p(q) [Eq. (12)) ] also has 
the form (|2"9"|) with the replacements Yli=i ~ > Si=i and 



ipfm) = n ci 



rn + I 'm + T 



■|0>, 



(42) 
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and the elementary neutral excitations are now related 
to p3i (q), P32(q), P4i(q), and p 42 (q) [see Fig. [3J. 

Again, one can show that the above four density oper- 
ators give rise to four independent boson operators, i.e., 



&l(q) 


= a«P3i(q), 


&i(q) = 


a 9 Pi3(- 


-q), 


4(q) 


= ^ q p 32 (q), 


6a (q) = 




-q). 


4(q) 


= a g p 4 i(q), 


63 (q) = 


Oiqp\i{- 


■q). (43) 


&l(q) 


= «gP42(q), 


6 4 (q) s 


a q p2i{- 


-q), 



which satisfy the boson algebra (|26p . Indeed, the com- 
mutator of each density operator with its correspond- 
ing Hermitian conjugate is also given by Eq. (|33|) with 
7=1,2 and J = 3,4. The expansion (|34f is replaced by 

pj/(q) = (PFM| j0/7 (q)|PFM)+^ 7/ (q) 

= ^^(Si,i+Si t 2)Sq,o + Sp n (q), (44) 

while Eq. (|3^|) is preserved, but now (7, J) = (1, 2) and 
(4,3). 

The set of creation and annihilation boson operators 
(1431) implies that Eqs. ([33)) should be replaced by 

Pll(q) = N^ q , - £ e-(^ 2 / 4 -^ k / 2 fel(q + k)fe 4 (k), 
k,i=l,3 

p 22 (q) = N^ q , - J2 e-(^ 2 / 4 - iqAk / 2 fel(q + k)fe 4 (k), 

k,i=2,4 

P33(q) = ]T e-^ 2 / 4+l i Ak / 2 fo 4 t (q + k)6 l (k), 

k,i=l,2 

p 44 (q) = e- (Z9)2/4+lqAk/2 6 4 t (q + k)6 l (k). (45) 

k,i=3,4 

Again, the expression (|2"5|) for the electron density oper- 
ator is preserved, apart from the changes Yli=i ~ * X)i=i 
and N0(5 q! o — * 2N0^ q ,o- When compared with the results 
of Sec. lII Cl the boson representation of the z-components 
of the spin, pseudospin, and mixed spin-pseudospin den- 
sity operators are interchanged, i.e., 

^(q) = £/T(q,k)&l(q + k)Mk), (46) 
i,k 

P z (q) = - e-^) 2 / 4 

x^cos(qAk/2)6t(q + k)6 i (k), (47) 

i,k 

with 7^(q) = 5z(q) and PzSz(q), and 

/f(q,k) = -/ 4 s (q,k)=*e-^ 2 / 4 sin(qAk/2), 
/f(q,k) = -/ 3 S (q,k) = e -^ 2 / 4 cos(qAk/2), 

and 

A PS (q,k) = -/f S (q,k) = - e -^ 2 / 4 cos(qAk/2), 
/ 2 PS (q,k) = -/ 3 PS (q,k) = -* e -^) 2 / 4 sin(qAk/2). 




FIG. 4: (color online) Schematic representation of the honey- 
comb lattice. The triangular sublattices A and B are respec- 
tively represented by blue and red circles, a and b are the 
primitive vectors of the underline triangular sublattice A and 
d is the basis vector. 

We again refer the reader to the Appendix [5] for the 
boson representation of the operators Pi2(q), P2i(q), 
p 34 (q), and p 4 3(q)- 

The generalization of the bosonization method-^ for 
the case of electrons restricted to the lowest Landau level 
and in the presence of two discrete degrees of freedom 
is concluded. The next sections will be devoted to an 
application of the formalism. 



III. QUANTUM HALL FERROMAGNETISM IN 
GRAPHENE 

In this section, we apply the methodology developed 
above to study the QHE at v = ±1 and v — in 
graphene. We will follow the lines of Ref.[ll| and derive an 
effective boson model for the system. Our starting point 
is the continuous model for graphene recently proposed 
by Alicea and Fisher^ Before outlining the derivation of 
this model, we will briefly review some aspects of the 
Landau level spectrum in graphene. 



A. Preliminaries on graphene 

Graphene is a collection of carbon atoms, which are 
arranged in a two-dimensional honeycomb lattice, as it 
is illustrated in Fig. \^r^ The lattice structure is trian- 
gular with two atoms per unit cell located at the posi- 
tions (0,0) and d = ao(0, l/\/3)- The lattice spacing is 
ao = 2.46 A. It might also be seem as two interpen- 
etrating triangular sublattices A and B. The primitive 
vectors of the (A) triangular lattice are a = ao(l,0) and 



8 



b = ao(— 1/2, v3/2), and therefore the primitive vec- 
tors of the reciprocal lattice are a* = (27r/a )(l, l/v3) 
and b* = (27r/ao)(0, 2/V3). In this atomic arrangement, 
the carbon atoms are connected by strong covalent a- 
bonds, derived from the sp 2 hybridization of the atomic 
orbitals. The remaining p z orbitals (perpendicular to 
the plane) have a weak overlap and therefore they form a 
narrow band of 7r-orbitals, through which the Fermi level 
passes.— By describing the 7r-electrons within a tight- 
binding model 



H t = -t 



(»!, 



(48) 



where t w 2.7 eV is the nearest- neighbor hopping en- 
ergy and the operators a\ a and b\ create a spin a 
electron on site i of the sublattices A and B respec- 
tively, one can show that the single-particle electron en- 
ergy varies linearly with momentum (e 9 = ±7u^p|q|, with 
vf = ao\/3t/2 » 10 6 m/s) around the six corners of the 
(hexagonal) Brillouin zone, i.e., the band structure con- 
sists of six Dirac cones. Only two of them are inequiv- 
alent, and here we consider the ones around the points 
K = (27r/o )(2/3,0) and K' = -(2 7 r/a )(2/3 J 0). In the 
undoped case, there is only one 7r-electron per carbon 
atom, the Fermi level lies at the Dirac points and there- 
fore the system is semi-metallic. By using a gate- voltage, 
it is possible to modify the carriers, either p-type or n- 
type (doped case). 

The fact that the electronic structure of the system 
may be described by an effective massless (continuous) 
Dirac model has some important consequences. In par- 
ticular, when a perpendicular magnetic field is applied, 
a different Landau level structure emerges when com- 
pared to the Schrodingcr-like one observed in the two- 
dimensional electron gas in GaAs heterostructures. In- 
deed, one can show that the energy of the (Dirac) Landau 
levels are given by 



E n . a = t\e z + sign(n)^2hv 2 F B\n\/c, (49) 

respectively for a = ] and |, which are associated with 
the two-components spinor eigenvectors (n 7^ 0) 



\g> \ - — ( l nm > 

\^ n ,m,a, a =+) - ^ y sign(n) |n - 1 m) 



(50) 



I* 



n,m,o".a — — / 



For n — 0, we have 



l$0,m,oy 



|0m) 




1 / sign(n)|n — 1 m) 

V l nm ) 



and |$ .„ 





|0m) 



(51) 

Here, a = ± corresponds, respectively, to K and K' 
points, m is the guiding center quantum number, and 
lnm) are the Schrodinger's Landau level eigenvectors. 



Each spinor component is related to one of the tri- 
angular sublattices A and B. For n ^ 0, each eigen- 
vector |$n,m,0-,+) and yt n ,m,<T,—) has a weight (proba- 
bility) equally distributed between the two sublattices, 
while the lowest Landau level eigenvectors \&o, m ,ir,+) and 
|^o,m,o-,-) are respectively localized on sublattices A and 
B. The results (|i9" ]) -([B"T ]) show that the Dirac Landau 
levels are approximately four-fold degenerate due to the 
electronic spin and valley (a — ±) degrees of freedom. 

Let us concentrate on the integer quantum Hall states 
in the lowest Landau level (n = 0). Apart from the fact 
that the fermion field operator Vl/(r) is a two-component 
spinor and the momenta q are measured with respect to 
the K and K' points [we refer the reader for a detailed 
discussion in the Appendix [C] , the methodology devel- 
oped in the previous section can be used to study the 
QHE at v = — 1 and v = 0. In fact, the former, which 
corresponds to a quarter filled lowest Landau level, is 
associate with the spin-pseudospin polarized phase (Sec. 
Ill Bp . whereas the latter, characterize by a half filled low- 
est Landau level, is associate with cither the spin (Sec. 
Ill C|) or the pseudospin (Sec. Ill D[) phases. 



B. Alicea and Fisher's model 

The effective continuous model proposed by Alicea and 
Fisher to study the quantum Hall effect in graphenegoes 
beyond the tight-binding approximation [see Ref. |8| for 
details]. In addition to 7i t [Eq. (T4"8|) ]. it also includes the 
on-site electron-electron repulsion term 7ijj and the (long 
range) Coulomb interaction 7icoui, namely 



H — Ht + Hjj + Hcoulj 



where 



Hu = Uj2 



i(n 2 ) 2 - ■ S, 



and 



Hcoul 



2^ 



V(ri - Tj)fknj 



(52) 



(53) 



(54) 



Here, U is the on-site repulsion energy and V(r) = e 2 /er 
is the Coulomb potential, with an estimated dielectric 
constant e w 5 [the energy scales for graphene are listed 
in the Appendix [S] . The electron number operator is 

fn = 4t c *T +4d c iU S * = ( 1 / 2 )E ct ,a4<tO'<tac 4 a, is the 
spin operator, where a is a vector of Pauli matrices, and 
c l a =a la or K a depending whether i is on sublattice A 
orB. 

Starting from the Hamiltonian (|52[) . a continuous in- 
teracting theory was derived by expanding the fermion 
operators a\ and b\ around the two Dirac points K and 
K'. After adding a perpendicular magnetic field B = Bz 
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and projecting into the lowest Landau level, the model 
may be rewritten as 



H = TLs 



U4 



where 



n su4 = ~^2v(q)p(q)p(-q) 



(55) 



(56) 



is the SU(4) invariant part of the Hamiltonian, with 
v(q) = 2ire 2 /eq (the Fourier transform of the Coulomb 
potential in two-dimensions), and 2 ^ 

Hsb = -^Sz(q = 0)-4^Ui(q)P z (q)Pz(-q) 



-p(q)p(-q)+P z (q)P z (-q) 



Is(q)-S(-q)-iPzS(q)-P z S(-q) 



(57) 



contains terms that break the SU(4) symmetry. The 
parameter uo is related to the on-site repulsion energy 
[no = VSoqU/A] and vi(q) is the Fourier transform of 



ui(r) 



V(r + y/V3) - (1 - 5 r , ) V(r)] . (58) 



The model (1551) was analyzed in two distinct situations: 
(i) the quantum Hall ferromagnetic regime, which cor- 
responds to an ideal, completely clean sample, and (ii) 
the quantum Hall paramagnetic regime, where disorder 
effects are very strong (very dirty sample). Here, we will 
only focus on the quantum Hall ferromagnetic regime. 
The interplay between disorder and electron-electron in- 
teractions will be postponed for a future publication. 

In order to derive an effective boson model for the 
quantum Hall states at v — —1 and v — 0, we just 
need to substitute the respective boson representation 
of the electron density, the spin, pseudospin and mixed 
spin-pseudospin density operators into the Hamiltonian 
(|55|) and normal order the resulting expression. Although 
the expansion of the electron density operator is simi- 
lar for the three phases, each quantum Hall state should 
be treated separately because the expansions in terms 
of bosons of the spin/pseudospin density operators vary 
from phase to phase. 



1. Filling factor v = — 1 

We start by considering the QHE at v = — 1 . The state 
at v = +1 is related to it by particle- hole symmetry and 
will not be discussed here. 

SU(4) invariant terms - Let us firstly analyzed the 
SU(4) invariant part of the Hamiltonian (|55|) . Substi- 
tuting the boson representation of the electron density 
operator [Eq. (|29|) ] in TLsua [Eq. (|56p ] and normal order- 
ing the boson operators, apart from a constant related 



to the positive background, we arrive at the following 
interacting boson model 



H 



SUA 



n 



n_/B 

Hi , 



where the quadratic part is given by 



< = EE w A f (qMq) 



(59) 



(60) 



and the quartic term reads 

3 

ijj=l q.P,k 

(61) 

The effective boson model (1591 is the SU(4) counter- 
part of the boson model derived in Ref. [Ill for the two- 
dimensional electron gas at v = 1 realized in GaAs het- 
erostructures (hereafter called 2DEG at v = 1). The 
ground state of the model ([59")) is the boson vacuum, 
which is the spin-pseudospin polarized state |SPFM). 
Notice that this state is indeed a spin polarized charge 
density wave (CDW) because the electronic distribution 
is concentrated only in one sublattice [see Eqs. (|5 1 1) and 
the discussion below this equation]. Hq describes three 
well-defined branches of bosonic excitations, character- 
ized by the same dispersion relation 



e~ 



(62) 



where Io(x) is the modified Bessel function of the first 
kindJ£ Eq. (l6"2")) is equal to the dispersion relation of 
the elementary neutral excitations (magnetic excitons) 
of the 2DEG at v = L n i 24 In the long wavelength limit, 
w q w £s|^q| 2 with e# = (e 2 /el), and therefore 

the branch i = 1 corresponds to spin wave excitations, 
while the branches i = 2 and i = 3 to pseudospin wave 
and mixed spin-pseudospin wave excitations, respectively 
[see Fig. [I]. At short wavelengths w q ~ \Zk/2 (e 2 /el), 
which is the energy of a very-well separated particle-hole 
pair. 24 Finally, the boson-boson interaction potential is 



given by 

u q (k,p) = 2v(q)e 



-(lq) 2 /2 



sin (qA k/2)sin (qAp/2) 



(63) 

Apart from the fact that Eq. (|6~lj) describes scattering 
processes between bosons within the same (i = j) and 
different (i ^ j) branches, the interaction potential (|6"3"|) 
is similar to the one derived in Ref. flTL It is worth men- 
tioning that our approach also provides an interaction 
between the bosonic excitations, which is not captured 
by the analysis presented in Ref. 0. 

Due to the similarities between the quantum Hall sys- 
tem in graphene at v = — 1 and the 2DEG at v = 1, we 
would expect that the charged excitations of the former 
might be described by topological soliton c 25 ! 26 (quantum 
Hall skyrmions) as well. In fact, the situation here is for- 
mally identical to the one in the (spinfull) bilayer QHS at 
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vt = 1 in GaAs heterostructures. The similarity clearly 
appears when the upper-layer and down-layer electronic 
states are combined into the bounding and anti-bounding 
states. In this case, there are four possible kinds of 
charged excitations with topological charge Qt = ±1 
(and corresponding electric charge Q e — eQr), namely 
one skyrmion (Qt = 1 and Q e — e) and three types 
of antiskyrmions (Qt — — 1 and Q e — —e). Indeed, 
they might be considered as SU(4) skyrmions because 
the topological excitation created by introducing an ex- 
tra electric charge should involve the three branches of 
neutral excitations in order to minimize the total energy 
[see Ref . I20I for a detailed description of SU(4) skyrmions 
in the context of the bilayers]. 

For the SU(2) version of the model we know that 
the boson-boson interaction potential (f63|) gives rise to 
bound states of two-bosons which are related to small 
[SU(2)] skyrmion-antiskyrmion pair excitations^ The 
fact that (|6"5|) describes scattering processes between dif- 
ferent bosonic branches indicates that here we would ex- 
pect bound states constituted by bosons belonging to the 
same and distinct branches. Moreover, it was also shown 
that by describing the topological excitation as a coher- 
ent state of bosons |sk) [see Eq. (64) of Ref. [ll|, the 
expectation value of the SU(2) boson model with respect 
to the state |sk) is equal to the energy functional derived 
from the phenomenological theory of Sondhi et al£^ for 
the quantum Hall skyrmion, i.e., the semiclassical limit 
of the SU(2) boson model agrees with Sondhi's theory for 
the quantum Hall skyrmion. 

It is quite straightforward to derive the semiclassical 
limit of the interacting boson model ([55]) . We start by 
writing down the SU(4) counterpart of the state |sk), 

|sk) =exp ^^^^6l(-q) + n^(q)^ |SPFM), 

(64) 

where (f2 l _ q )* = il q and the constant B will be de- 
termined latter. With the aid of the Baker-Hausdorff 
formula, one can show that the expectation value in 
the state |sk) of normal ordered boson operators is ob- 
tained just by replacing each &J(q) and &,(q) respec- 
tively for i23f2 q and — iSfi!_ q . Defining the excess charge 
Sp(q) — (sk\p(q)\sk) — N^t^o, where p(q) is the electron 
density operator (j!?9"j) . we have 

Sp(q) = 2z ]T £ er{lq)2/2 sin (1 A k /2) ^k +q f4- (65) 

i=l k 

Assuming that the Fourier transform of f2 q and f2 q 
vary slowly in space, we can restrict ourselves to the 
long wavelength limit of Eq. (|65p . i- e -j we can consider 
e -(k) /2 s j n (q /\ k/2) w qAk/2. Within this approxima- 
tion, the Fourier transform of 6p(q) is given by 

Sp(r) = iB 2 l 2 z ■ Vfi*(r) x Vfi l (r), (66) 



which is in agreement with the expression for the topo- 
logical charge density derived by Arovas et al£Z- in their 
studies of SU(N) quantum Hall skyrmions. Indeed, by 
comparing Eq. (f6"6"|) with Eq. (3) from Ref. H3, one con- 
cludes that B = I/V^ttP. Once the constant B is fixed, 
we can now calculate (sk\Ti.g U4 \sk) and show that 

(Wfw) = 2p° 5 ]T / d 2 r|V^(r)| 2 

+ \j d 2 rd 2 r'v(\r-r'\)Sp(r)Sp(r r ), (67) 

where = l/(16\/27r) e 2 /el is the stiffness and v(r) = 
e 2 /er is the Coulomb potential. The energy functional 
^[^(r)] = (Hg U4 ) [Eq. ([B7])]. which corresponds to the 
SU(4) counterpart of Sondhi's model, agrees with the 
findings of Arovas and co-workers 31 This analysis shows 
that the boson model ([59]) can indeed be used to study 
SU(4) quantum Hall skyrmions in graphene. 

Symmetry breaking terms - The degeneracy of the three 
branches of boson excitations is lifted when the SU(4) 
symmetry breaking term 7is,g is taken into account. Fol- 
lowing the same procedure used above, we can derive an 
effective boson model from the Hamiltonian (|57|) . The 
task here is slightly more difficult because TLsb involves 
more complex expressions. 

Let us start by expanding the operators S(q) and 
PzS(q) in terms of the density operators pjj(q). It is 
possible to show that 

S(q) ■ S(-q) + P z S(q) • P Z S(-q) = 
S z (q)Sz(-q) + PzSz(q)PzSz(-q) + Pi2(q)p2i(-q) 

+ P2i(q)pi2(-q) + P34(q)p43(-q) + P43(q)p34(-q)- 

The boson representations of the above density operators 
Pu(q) are shown in the Appendix IB1 [see Eqs. (|B1[) and 
(|B2[) ]. After a lengthy but straightforward calculation, 
one can show that 7isB is also mapped into an interacting 
boson model. Adding Eq. (|59p. which was derived from 
the SU(4) invariant term, the total effective boson model 
may be written as 

H B =n^+Hf. (68) 
The quadratic term now reads 

i— 1 q 

where u>i(q) are the renormalized boson dispersion rela- 
tions, 

= Pz + 2( Uo - U i)N (l-e-^ 2 / 2 ) +w q , 
uj 2 (q) = 4u 1 N -4^w 1 (k)e- ( " c)2/2 

k 

x cos 2 (k A q/2) + w q , 
w 3 (g) = E z + oj 2 (q), (70) 
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with w q given by Eq. (f62|) and Ux — ui(q = 0) [see Eq. 
(|58[) ] . In the small momentum region, we have 

wife)* A i + 47r4|Zq| 2 , (71) 

where the excitation gaps Aj and the renormalized stiff- 
nesses are given by 

Ai = E z , A 2 = V7r 3 /24(ao/OwiN , 

A 3 = A! + A 2 , 

= N^(«o-ui)/47r + pS. 
p| = p| = Ul N /4^ + p^. (72) 

Notice that ui2(q = 0) <wi(g = 0) and uj^q = 0). Both 
small and large momentum limits of Eqs. (|70[) agree with 
the results derived by Alicea and Fisher.— The interaction 
term assumes the form 

3 

= E E p )^ k + q ) 6 K p - q)^(p)^( k )' 

(73) 

where the total boson-boson interaction potential 
■Dq J (k, p) , which is richer than the one derived only from 
the SU(4) invariant part of the Hamiltonian is given 
by 

tJ^(k,p) = ^-^e^/ 2 (sin (q A k/2) sin (q A p/2) 
- cos (q A k/2) cos (q A p/2)) 

+ ^,1^.1 ^e-' 2 li +k l 2 / 2 cos (q A (k - p)/2) 

+ 4 (^-u 1 )/f(q,k)tf(-q,p) 

+ kj ^-e- (lq)2/2 (2 sin (q A k/2) sin (q A p/2) 

+ e iqA(k-p)/2^ 

+ ^£ e - q A(k- P )/ 2 r 6 . i(1 _ Me ^| q+k p/2 

3 V 

- 4 2 ^ 3 e-' 2|q+k - p|2/2 ) , (74) 

with the form factors ff (q, k) given by Eqs. (J2J) and 
5ij = 1 — 5ij. Finally, it is worth mentioning that the 
ground state of the system is still the boson vacuum 
|SPFM). Indeed, this result is corroborated by exact 
diagonalizations on small systems^ 

The introduction of new terms in the boson-boson in- 
teraction potential might modify the two-bosons spec- 
trum, for instance, one particular kind of bound state 
may have lower energy than the others. As a conse- 
quence, one specific type of skyrmion-antiskyrmion pair 



excitation will be more favorable. Indeed, it was ar- 
gued that a pseudospin skyrmion-antiskyrmion excita- 
tion should determine the charge gap due to the small- 
ness of the excitation gap tl>2 (5 = 0)£ 

Although v q J (k, p) is quite complex, it is possible 
to make a simple analysis by considering the SU(2) 
limit of the bosonic Hamiltonian (|68|) and then by 
calculating the semiclassical limit of this reduced bo- 
son model. In this case, Eq. (|64|) simplifies to 
|sk) = exp(B£ q Q q &J(-q) + n q 6 4 (q))|SPFM), i.e., we 
assume that only the i — th bosonic branch is exci- 
tated while the others are kept frozen. Following the 
same steps which lead to Eq. (|67[) and approximating 
e -0?) 2 /2 cos ( qA i c /2)cos(qAp/2) « er^ 2 ' 2 , one can 
show that the functional energy for the different skyrmion 
flavors assumes the form 

E % [n(r)] = [n] + Ef [n] + Ef z [n] + Ef [n] 

= J d 2 r[2£ 2 (l-A 4 n z (r))+2/^(Vn(r)) 2 ] 

x (l-n z (r))(l-n z (r')) 

+ \S '< iW (^+^-"-'' l " /2 '") 

(75) 

where i — 1,2,3 refer respectively to spin, pseudospin, 
and mixed spin-pseudospin-like skyrmions. n(r) is a unit 
vector defined by the relation f2(r) = z x n(r) [see Ref. 
|Tl| for details]. Aj and p l s are given by Eqs. (j?2"j) . fij = 0, 

u 2 = u 3 = 32N0(u o /6 - Ui), v\ = 2N0(2u o /3 - Hi), and 

V2 = V3 = N0M o /3. 

Notice that the SU(4) symmetry breaking part of the 
Hamiltonian (|55f adds to the energy functional for the 
skyrmion a Zeeman-like term (i^ A [n]) and provides small 
contributions to both the stiffness and the topological- 
charge-topological-charge interaction potential. In fact, 
V2 and V3 > while v% can be either positive or negative 
depending on the value of the on-site repulsion energy 
U. For the pseudospin and the mixed spin-pseudospin- 
like skyrmions (boson branches i = 2 and 3, respec- 
tively) there is an extra contribution given by Ef z [n]. 
This term favors excitations with n z (r) < because 
4N0(m o /6 - mi) < when 2 < U < 12 eV [see Appendix 
[A] . Remembering that a quantum Hall skyrmion is char- 
acterized by n z (r) — > 1 when r — > oo, one might conclude 
that E zz [n] contributes to an increase of the radius of 
the skyrmion and so its stability. 

It is still difficult to predict which type of skyrmion 
has the lowest energy without performing careful calcu- 
lations. What we can easily see is that if the on-site re- 
pulsion energy is U ss 10 eV, then the scenario proposed 
in Ref. [H, that a pseudospin skyrmion should be the low- 
est energy one, is confirmed. In this case pg ~ p 2 s ~ p s 
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and v\ ~ V2 ~ £3 <C e c . As the skyrmion energy is 
47rpg + C(A i /ec , ) ) 25 i 26 the lowest energy soliton should 
be the one related with the excitation branch which has 
the smallest Aj, i.e., the pseudospin branch (i = 2). No- 
tice that the presence of Ef z [n] for i — 2 and 3 does 
not alter the above conclusions because this term should 
reduce the total energy. 

Finally, it is worth mentioning that (|74p contains 
cos (q A k/2) cos (q A p/2) like terms, which are also 
present in the boson-boson interaction potential derived 
for the bilayer QHS at Vt = 1 (spinless case) within 
the SU(2) bosonization method. 15 Such similarity im- 
plies that, in principle, a Bose-Einstein condensate could 
be realized here. Let us consider again the SU(2) limit 
of the boson model |68|) and focus, for instance, on the 
mixed spin-pseudospin branch (i = 3). The phase with 
N^/2 bosons should then correspond to the antiferro- 
magnetic one proposed by Herbut^ Assuming that the 
bosons condense in their lowest energy mode (q = 0) 
and treating the reduced boson model within the Bo- 
goliubov approximation, one arrives at a model simi- 
lar to Eq. (8) from Ref. [TH with the replacement 
A q -»■ 8(m /6 - m) exp(-(7q) 2 /2). As it was showed in 
the last paragraph, uq/6 — U\ < 0, and therefore such a 
phase should be unstable. Similar considerations hold for 
the pseudospin branch. The situation is more delicate for 
the spin wave branch and it will not be discussed here. 



2. Filling factor v — 

The analysis of the quantum Hall state at v — follows 
the same lines of the previous section with the difference 
that now either the spin or the pseudospin phases can be 
realized. 

Spin phase - Let us firstly assume that the system 
is in the spin phase. In this case, an effective boson 
model can be obtained from the fermionic Hamiltonian 
([55)1 with the aid of the expressions calculated in Sec. 
ILTCl and Eqs. (JB3J) and (|B"4]) . Due to the fact that 
the boson representation of the electron density opera- 
tor (fi"2")) does not change from phase to phase, the boson 
model derived from the SU(4) invariant part of the to- 
tal Hamiltonian [Eq. (|56p] is similar to Eq. ([55]) with 

the replacement Yl%=i ~* 2<=i* The ground state of 
the system is also the boson vacuum, which now corre- 
sponds to the state |SFM). There are four branches of 
well-defined bosonic excitations. In the small momentum 
region, the branches i = 1 and 4 describe spin wave exci- 
tations whereas the branches i — 2 and 3 correspond to 
mixed spin-pseudospin wave excitations [see Fig. [2]. 

The total effective boson model, which also includes 
the terms obtained from TLs b , reads 



H B = U B + nf + V B . 



(76) 



The quadratic term H B is again given by Eq. (f6"5)) with 

wi(g) = u) 4 (q) 



= E z + 2(u - Ul )N* (l - e -(^) 2 /2) + ^(77) 
u> 2 (q) = u 3 (q) 

= E z + 2w o N - 2 Ul N (l + e -<^ 2 / 2 ) + w q , 

In the small momentum 



where w q is given by Eq. 
region, Qi(q) assume the form (|7ip . with the following 
excitation gaps A, and stiffness p l s 



Ai = A 4 =E Z , 

A 2 = A 3 = E z + 2N,i(uo - 2«i), 

Ps = Ps = N ^( u o -ui)/4tt +p° s , 

Ps = Ps = N ui/4tt + p^. 



(78) 



Notice that the introduction of the symmetry breaking 
terms does not modify the ground state of the system 
|SFM). 

The boson-boson interaction part of the total Hamilto- 
nian has two distinct terms. The first one, ftf , is equal 
to Eq. (|T3"|) , but now the interaction potential v 1 ^ (k, p) 
reads 



<(k, P ) 



u e 
1 



{lq)2/2 (sin (q A k/2) sin (q A p/2) 
cos (q A k/2) cos (q A p/2) 



" kjf (e- ( ' 9)2/2 - 2(^,1 + M 

xe- ;2|q+k|2/2 ) cos(qA (k-p)/2) 
+ 4(|- Ul )/f(q,k)/f(-q )P ) 
+ hi ? \hi ,]Z~ (lq)2 ' 2 cos (q A (k - p) /2) 
-^' 2 sin (qA (k-p)/2) 



e -i 2 |q+k| 2 /2 e - 4 (-l) l +^qA(k-p)/2 



, (79) 



with the form factors ff (q, k) given by Eqs. (|4Tj) . hij = 

^i,i«5j,4 + <*»,4<5i,i + <5i,2^,3 + ^i,3^,2 ; and hij = 1 - hij 
The second one, V s , can be written as 

VB = E <(k,p)&I(k + q)&i(p-q)fe 3 (p)6 2 (k), (80) 
q,p,k 

where 



<(k,p) = ^%xp(zqA(p-k)/2) 



x (gC-^k+q)! 2 +e -|i(p-q)l" 



(81) 



Pseudospin phase- Turning to the pseudospin phase, 
similar considerations show that this phase is also char- 
acterized by an effective boson model analogous to (|T6[) . 
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The ground state is the boson vacuum |PFM) [Eq. (|42f ] 
and the dispersion relations of the four branches of 
bosonic excitations are 



uji(q) = E z (5 ir 3 - 5 i)2 ) - 2u N,s, 
+ 2 Ul N (3 - e -('?) 2 /2^ + . 



(82) 



for i = 1, 2, 3 and 4. The long wavelength limit behavior 
of uii (q) is also given by Eq. (fTTj) with 

Ai = A 4 = 2N (2ui-m o ), 

A 2 = -E z + A u A 3 = £' z +A l! (83) 
^ = N^ui/47r + p&, * = 1,2,3,4. 

Here the branches i = 1 and 4 describe pseudospin wave 
excitations, while i = 2 and 3 correspond to mixed spin- 
pseudospin wave excitations [see Fig. [3J. The boson- 
boson interaction potential in Tif is 



i#(k,p) 



(uq sin (q A k/2) sin (q A p/2) 



(uq - 4ui) cos (q A k/2) cos (q A p/2)) 
^.!^e-^ 2 / 2 cos(qA(k-p)/2) 



3 

U 



J %3 



^,ie~ ( ' 9)2/2 cos(qA (k-p)/2) 
_iV'+i P -('9)72, 



/Wj (i(-l) l+J e- wr/ " sin (q A (k - p)/2) 



e -; 2 |q+k-p| 2 /2 e -i(-l) i+3 qA(k-p)/2 



(84) 



Finally, the interaction term V B can be written as 

vB = E "%(k,p)(&t(k + q)6|(p-q)6 3 (p)b 2 (k) 
q,p,k 

+ 4(k + q)^(p-q)& 4 (pMk)) , (85) 



with 



xe 



2wo 
3 

-iqA(p— k)/2 



„» q (k, P ) = _^£( e -^) 2 /2 +e -Wq+k-p|) 2 A 

3 V 



(86) 



The small and large momentum expansions of Eqs. 
(ITT)) and l|82p are in agreement with the results of Al- 
icea and Fisher^ who present a detailed discussion about 
the stability of each phase. Here we just want to point 
out that the behavior of the smallest excitation gap in- 
dicates which phase should set in. For instance, in the 
spin phase, u)2{q = 0) and 0)3(9 — 0) are smaller than 
u)\{q = 0) and u>±{q = 0) as long as uq — 2u\ < 0. This 
result implies that the spin phase is stable only if 

< a> 2 , 3 (g = 0) = E z + 2N (w o - 2ui), 



where the estimated values of the parameters Ez, uq, and 
u\ are shown in the Appendix fAJ It is possible to show 
that the spin phase sets in only if U > Uc ~ 3.25 eV. 
The opposite condition is found by carrying out the same 
analysis in the pseudospin phase. It is difficult to con- 
clude which phase is more favorable due to the uncertain- 
ties in the determination of the on-site repulsion term U. 

Charged excitations - Concerning the elementary 
charged excitations, the similarities between the effec- 
tive boson model derived from Hsua [Eq. (|56")) ] and the 
SU(2) counterpart^ indicate that, in both phases, the 
lowest energy charged excitations should be described by 
quantum Hall skyrmions as well. Again, within our for- 
malism, such kind of topological excitation is given by 
the state (|6"4")l . This scenario agrees with the numerical 
calculations of Yang et al.^ who show that skyrmions 
should occur in the n = as well as n = 1, 2, and 3 
Dirac Landau levels. 

We can carry out the same analysis of the previous 
section and calculate the semiclassical limit of the cor- 
responding SU(2) boson models for the spin and pseu- 
dospin phases in order to estimate how Hsb changes 
the skyrmion energy. It is easy to show that the en- 
ergy functional i?i[ n ( r )] is as in Eq. (|75|) . For the spin 
phase, the parameters of -E;[ n ( r )] are u\ = U4 = 0, 
U2 = u 3 = 32N0(u o /6 - ui), V\ = Vi = 2N <j) (2u /'i - u\), 
v-2 = V3 = N0ito/3, and Aj and p l s are given by Eqs. (|78[) . 
Again, if U » 10 eV, the corrections due to the SU(4) 
symmetry breaking terms are such that all stiffnesses are 
equal and so the topological-charge-topological-charge in- 
teraction potential. In this case, we also have Ai = A4 w 
A2 = A3, and therefore both mixed spin-pseudospin and 
spin skyrmions can be realized. 

For the pseudospin phase, £?j[n(r)] is characterized 
by mi = u 2 = u 3 = in = 32N0(w o /6 - «i = 
v-2 = V3 = V4 = N0Mo/3, and Ai and p l s given by 
Eqs. (|55| . Here the SU(4) symmetry breaking terms 
equally affect the four excitation branches, independently 
of the value of the on-site energy U . The mixed spin- 
pseudospin skyrmions should have a larger radius (and 
therefore lower energy) than the pseudospin ones because 
A2/e c < Ai/e c and l±i/e c as it was already pointed out 
in Ref. [8j. The term E zz [n(r)] will even reduce the en- 
ergy of the texture with larger radius. 

Finally, we should emphasize that our bosonization 
method is designed to study only bulk excitations. The 
recent proposal of Abanin and co-workers^ that the fi- 
nite value of the longitudinal conductivity at v = is 
related to the existence of charged gapless excitations at 
the edge of the system can not be addressed with our 
formalism. 



IV. SUMMARY 

We presented here a non-perturbative bosonization 
scheme for electrons restricted to the lowest Landau level 
in the presence of two discrete degrees of freedom, spin- 
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1/2 and pseudospin-1/2. We analyzed the cases when the 
lowest Landau level is quarter-filled and half-filled. For 
the latter, two distinct phases can be realized, the so- 
called spin and pseudospin phases whereas in the former 
only the spin-pseudospin phases sets in. In each case, a 
set of n-independent kinds of creation and annihilation 
boson operators were defined and the boson represen- 
tation of the projected electron, spin, pseudospin, and 
mixed spin-pseudospin density operators were calculated. 
The bosonic expressions derived obey the lowest Landau 
level algebra. 

We then applied the formalism to study the QHE at 
v = and v = —1 in graphene. We concentrated on very 
clean samples, assuming that the system is in the quan- 
tum Hall ferromagnetic regime. For each quantum Hall 
state, the continuous fermionic model proposed by Alicea 
and Fisher 8 was mapped into an effective interacting bo- 
son model. We showed that the quadratic term of this 
model describes n well-defined branches of bosonic exci- 
tations, whose dispersion relations are in agreement with 
the asymptotic ones calculated by Alicea and Fisher. 8 
Our formalism allows us to go beyond the analysis pre- 
sented in Ref. |8| as we are able to calculate the interaction 
between the n bosonic excitation branches. 

The boson model Tig U4 derived from the SU(4) in- 
variant part of the fermionic Hamiltonian is similar to 
its SU(2) counterpart obtained before in our studies of 
the 2DEG at v = 1. Based on this analogy, we argued 
that the charged excitations for the quantum Hall states 
in graphene should be describe by topological solitons 
(quantum Hall skyrmions) and proposed that such exci- 
tation can be written as a bosonic coherent state |sk), 
generalizing the SU(2) expression of Ref. [TJ We then 
calculated the semiclassical limit of Hgjj^ and showed 
that the derived energy functional is equal to the one 
calculated by Arovas et al. for SU(N) quantum Hall 
skyrmions j22 

We briefly discussed how the SU(4) symmetry breaking 
terms modify the skyrmion energy functional by taking 
SU(2) limits of the total boson model and then calcu- 
lating the semiclassical limit of the reduced models, i.e., 
focusing on one specific skyrmion flavor. We showed that 
both the stiffness and the topological-charge-topological- 
charge interaction potential are renormalized and that an 
extra term (Ef z [n\), which favors larger skyrmions ra- 
dius, is introduced. More detailed studies of the boson- 
boson interaction potential as well as of the disorder ef- 
fects are deferred to a later publication. 

The method presented here is quite general. It might 
be used to study bilayer quantum Hall systems at vt = 1 
and vt — 2 realized in GaAs heterostructures. In partic- 
ular, it will allows us to address questions related to the 
electronic spin, which seems to play an important role in 
the behavior of these systems at vt = L 
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APPENDIX A: ENERGY SCALES 

The relevant energy scales for the QHE in graphene 
are presented in Table [A"l The cyclotron (hwc), Zee- 
man (Ez), and Coulomb energies (ec) as well as the 
parameters N^uo and N^iii [see Sec. IIII B| are given 
in terms of the magnetic field B, which is measured in 
Tesla. We consider the following estimated parameters 
for graphene: effective g-factor g ss 2, dielectric constant 
e « 5, and on-site repul sion ener gy 2 < U < 12 eV The 
magnetic length I = \fhcj (eB) = 256 /yB is measured 
in angstroms and ao = 2.46 A is the lattice spacing of the 
triangular underlining lattice. 

TABLE I: Energy scales for the QHE in graphene 



Energy scales 


(K) 


hwc 


y/2hv%B/c 


380.60^ 


E z 




1.0&B 


ec 


e 2 /d 


150.1271? 


N^uo 


V3Ual/8Tvl 2 


0.08(7B 




aolec/V3 


0.4B 



APPENDIX B: BOSON REPRESENTATION OF 
THE DENSITY OPERATORS pi 2 (q) p2i(q), P34(q), 
AND p 43 (q) 

Let us concentrate on the spin-pseudospin phase. Al- 
though the boson operators &i(q) and b\(q) are respec- 
tively defined by pi2(q) and /02i(q), the boson repre- 
sentations of these density operators are not necessarily 
pi 2 (q) = Q!q 1 6i(-q) and p 2 i(q) = a^b'^q). If it were 
the case, we would have [pi2(q), P2i(k)] = ar <5 q ,_k, in 
completely disagreement with the commutator (|SJ). The 
same procedure described in Sec. IIIBI should be em- 
ployed in this case as well. Due to the similarities be- 
tween the steps involved here and in the calculation of 
the boson representation of the spin density operators 
S" + (q) and 5 ,_ (q) of the SU(2) case, we refer the reader 
to Sec. II. C of Ref. [U for all the details and just display 
the final results here. We have, 

p 2 i(q) = ^e-^ 2 / 4 &t(q), 

(Bl) 



15 



- J2 // 2 (q,k,p)6 4 t (k + q + p)&. 1 (p)6 1 (k), 

i,k,p 

where the form factors are given by 

A 12 (q,k,p) = N <A - 1 /2 e -fe) 2 /4 cos((q + k)A(p + q)/2)) 
/ 2 12 (q,k,p) = / 3 12 (q,k,p) 

= N _1/2 e~ (i9)2/4 e~ j(q+k)A(p+q)/2 , 

and 

P3i(q) = P34(-q) 

= e -('9) 2 /4^ e .qAk/2 6 t (q + k)&3(k ^ (B2) 
k 

Similar considerations hold for the spin phase. In this 
case, the boson representation of both P2i(q) and P43(q) 
arc defined respectively by the creation boson operators 
b\(q) and b\(q), while more involve expressions are de- 
rived for their Hermitian conjugates. We have, 



p 21 (q) ee VN^-^/V^q), 



(B3) 



- ]T/ 34 (q,k,p)^(k + q + p)6 2 (p)& 3 (k), 
k, P 

with the following form factors 
/ 2 34 (q,k,p) = / 34 (q,k,p) 

= N ~ 1/2 e~ (i9)2/4 e~ j(q+k)A(p+q)/2 , 
/ 3 34 (q,k,p) = N0 -i/2 e -(i< Z ) 2 /4 e + l (q+k)A( P+q )/2 j 

/ 34 (q,k,p) = N^-^e-P'^oos^q + kjAtp + q)^). 

Finally, the pseudospin phase is characterized by the 
following expressions 

p 12 (q) = -e-^/^e-^/^ft^q + ^bxCk) 



6l(q + k)6 3 (k) 



PsM) = e-(^ 2 / 4 ^e iqAk / 2 (6 t 1 (q + k)6 3 (k) 

k 

+ 6 2 (q + k)fe 4 (k)), 



(B5) 



with p 2 i(q) = /f* 2 (-q) and P43(q) = peUC-q)- Thc 

derivation of the above expressions is analogous to the 
one involved in the calculations of Eq. (|B2[) . 



- E E z' 2 ^ k ' p) & l( k + «i + p)^(p) & i( k ) 

i— 1 k,p 

- - ^ / 12 (q, k, p)ftt( k + q + p)& 3 (p)6 2 (k), 

k,p 

where the form factors are 

A 12 (q,k,p) = N <A - 1 / 2 e-^ 2 /4 cos((q + k)A ( p + q)/2); 
/ 2 12 (q,k,p) = / 12 (q,k,p) 

= N«,4 _:L/2 e~ (i9)2/4 e +l(q+k)A(p+q)/2 , 
/ 3 12 (q,k,p) = N ^i/2 e -(i9) 2 /4 e -,(q+k)A(p +q )/2 ; 
and 

p 43 (q) ee ^e^ 2 /%(<D, 



(B4) 



P34 



(q) = V^^ )2 /%(-q) 



E fi 4 (^ k ' p) fe l( k + q + p)MpMk) 

i=2 k,p 



APPENDIX C: BOSONIZATION AND DIRAC 
LANDAU LEVELS 

For electrons in graphene subject to a perpendicu- 
lar magnetic field, the fermion field operator is a two- 
component spinor, which may be written in Dirac Lan- 
dau level basis as 



(CI) 



where K = (27r/ao)(2/3, 0) and the spinors |i> n ma) are 
given by Eqs. (|50|) and (fSTjl . 

Defining thc density operator paa,/3x(^) as 



(C2) 



one can calculate its Fourier transform in the same way 
as it is done in Eq. (J5]), i.e., 

p a *.M<D = E E^«-«i e " i(q+(a ^ )K) ' r i^'-'^) 

n,n' m,m' 

XC nma<7 C n' m' (3 A ■ (C3) 
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The projection into the n-th Dirac Landau level is ob- 
tained by taking the component n — n' in Eq. (|CJ3[> . For 
the lowest Landau level, the fact that the eigenvectors 
(|5ip have only one non-zero entry implies that 



-(lq) 2 /2 F a{3 
r n=0 



(C4) 



where F"f = S a ,p and /W„3A(q) is given by Eq. (5]). In 
the isospin language, p/j(q) = /0/j(q) and does not van- 
ish only if (I, J) = (1,1), (1,2), (2,1), (3,4), and (4,3). 
Notice that these are the density operators which appear 
in the effective continuous model (|55[) . Therefore, the 
expressions derived in Sees. Ill Bl - III Dl can be directly 
employed to study the fermionic model (|55|) . 

The situation is more complex for higher Landau levels 
(n 7^ 0). From Eqs. (|50|) . it is possible to show that^S 



puir,f3\{q) 



iq-!K| 2 /2 F a/3( q ) 



(C5) 



where the form factors F^^(q) read 



^? Q (q) 



L\ n \-1 



({kf 



\ 2 



(C6) 



(lq x - IK X ) Ll 



1-1 



IK\ 



and F„f + (q) = (— q) with K — * K. In the expres- 
sions above, we used the fact that 2K = K' = — K. The 
connection with the formulae derived in Sees. IIIBI - HTPl 
is obtained via the relation. 



Paa,j3\ 



(Zq - IK) = F 7 f(lq - ZK)p Q(Ti/9A (Zq) 



(C7) 



We refer the reader to Ref. |3C| for a detailed analysis 
of the form factors F" /3 (q) and theirs implications in the 
dynamics of the quantum Hall states in higher Dirac Lan- 
dau levels. 
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